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ABSTRACT 


Mathematical  tools  were  developed  previously  to  analyze 
the  dynamic  behavior  of  tethered  balloon  systems.  The  model  for 
the  tethered  balloon  system  consists  of  the  streamlined  balloon 
and  a  tether  made  up  of  three  discrete  links.  The  tools  consist 
of  the  linearized  characteristic  equations  which  incorporate  the 
physical,  aerodynamic  and  mass  characteristics  of  the  system  and 
the  dynamic  simulation  computer  program  which  determines  the 
response  of  the  tethered  balloon  system  to  wind  disturbances.  A 
tethered  balloon  system  consisting  of  a  70,000  cubic  foot  aero- 
dynamically  shaped  balloon  and  a  0.52  inch  diameter  Nolaro 
tether  was  flown  at  Fair  Site  on  White  Sands  Missile  Range  to 
obtain  experimental  motion  data.  Experimental  test  conditions 
and  test  tethered  balloon  system  characteristics  were  input  in 
the  computer  programs  and  balloon  dynamic  behavior  was  predicted. 
A  comparison  of  experimental  and  predicted  dynamic  characteris¬ 
tics  of  the  tethered  balloon  showed  a  reasonable  match  although 
some  discrepancies  existed.  In  spite  of  these,  it  is  concluded 
that  the  stability  and  dynamic  simulation  computer  programs  are 
a  powerful  tool  for  design  and  analysis  of  tethered  balloon 
systems . 
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SECTION  I 


INTRODUCTION 


The  objective  of  the  present  program  is  to  compare 
actual  flight  data  with  analytically  predicted  data  for  a 
tethered  balloon  and  to  improve  the  mathematical  model  if 
required.  A  second  objective  is  to  study  the  dynamic  behavior 
of  two  tethered  balloon  types 

Mathematical  tools  have  been  developed  on  a  previous 
program  to  analyze  the  dynamic  behavior  of  tethered  balloon 
systems.  The  techniques  used  are  determination  of  the  roots 
of  the  linearized  characteristics  equations  which  incorporate 
the  physical,  aerodynamic,  and  mass  characteristics  of  the 
system,  and  dynamic  simulation  of  the  tethered  balloon  systems 
to  determine  response  of  the  systems  to  wind  disturbances. 

The  techniques  are  complementary  and  each  helps  to  obtain 
insight  into  the  behavior  of  tethered  balloon  systems. 

The  model  for  the  tethered  balloon  system  consists  of 
the  streamlined  ballooii  and  a  tether  made  up  of  three  dis¬ 
crete  links.  The  derivation  of  non-linear  equations  of  motion 
for  this  system  were  devised  in  three  dimensions.  The 
equations  are  linearized  for  stability  analysis  and  treated 
as  uncoupled  in  the  longitudinal  and  lateral  degrees  of  free¬ 
dom.  Characteristic  equations  of  the  system  are  developed  and 
solved  for  the  roots  which  represent  the  frequency  and  damping 
qualities.  References  1,  2,  and  3  document  the  results  of 
this  program. 

A  tethered  balloon  system  consisting  of  a  70,000  cubic 
foot  aerodynamically  shaped  balloon  and  a  . 52  inch  diameter 
Nolaro  tether  was  flown  by  AFCRL  at  White  Sands  Missile  Range 
to  obtain  experimental  balloon  motion  data.  These  experimental 
data  have  been  compared  with  those  predicted  by  the  mathemat¬ 
ical  model.  The  efforts  to  establish  the  validity  of  the  anal¬ 
ytical  techniques  are  reported  herein. 


SECTION  II 


TECHNIQUES  FOR  STUDY  OF  DYNAMIC  BEHAVIOR  OF  TETHERED 

BALLOONS 


A.  MATHEMATICAL  MODELS 

A  system  of  differential  equations  was  developed 
(see  References  2  and  3)  that  describes  the  motion  of  the 
tethered  balloon  in  three  dimensions.  The  degrees  of  freedom 
associated  with  the  motion  are  yaw,  pitch  and  roll  of  the 
balloon  about  its  dynamic  mass  center,  and  pitch  and  yaw 
(lateral  rotation)  of  the  tether.  There  are  a  total  of 
3  +  2N  degrees  of  freedom  where  N  is  the  number  of  links  used 
to  simulate  the  tether. 


First  consider  the  longitudinal  degrees  of  freedom. 
The  dependent  variables  shown  in  Figure  1  are  9  (pitch  of  the 
balloon)  and  £  (pitch  of  the  "r"th  link) ,  where  r  is  a 
particular  lirR.  All  angles  are  shown  positive. 


Figure  1.  Balloon  Tether  Model  in  Longitudinal  Plane 


Preceding  page  blank 
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In  Figure  1  VfiR  is  the  relative  velocity  of  the 

balloon's  center  of  gravity  with  respect  to  the  air  and  is 
the  resultant  of  the  steady  wind,  the  wind  gust,  the  balloon 
translational  motion  and  the  velocity  due  to  rotation  of  the 
balloon  about  its  center  of  mass.  The  angle  of  attack  (a)  is 
the  angle  that  the  relative  wind  forms  with  the  longitudinal 
axis  of  the  balloon. 

The  lateral  degrees  of  freedom  are  displayed  in 
Figure  2  which  gives  the  front  and  top  view  of  the  tethered 
balloon.  The  lateral  degrees  of  freedom  are:  ifi(yaw  of 
balloon),  <f>  (roll  of  balloon),  and  o  (yaw  of  "r"th  link). 

All  angles  are  shown  positive) . 


Figure  2.  Balloon  Tether  Model  in  Lateral  Plane 


4 


Pertinent  geometry  of  the  tethered  balloon  and 
applied  forces  are  identified  in  Figure  3. 


In  order  to  separate  the  equations  of  motion  into  a 
longitudinal  response  and  a  lateral  response,  it  was  further 
assumed  that  the  system  was  near  equilibrium.  This  resulted 
in  a  set  of  equations  describing  the  longitudinal  motion 
which  is  coupled  only  in  the  pitching  variables  of  the  balloon 
and  the  pitching  variables  of  the  tether.  However,  the  second 
set  of  equations  for  the  lateral  motion  does  not  completely 
uncouple  from  the  longitudinal  degrees  of  freedom  because  the 
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-  Gross  static  lift  (Total  weight' 

of  displaced  air) 

-  Total  weight  of  balloon  including 

internal  gases 

-  Vertical  aerodynamic  forces 

acting  on  balloon 

-  Horizontal  aerodynamic  forces 

acting  on  balloon 

-  Aerodynamic  pitching  moment 

acting  on  balloon 

-  Aerodynamic  reference  center 


CB  -  Center  of  buoyancy  of  total  displaced  volume 

Cm  -  Apparent  mass  center  and  dynamic  center  in  pitch 

(and  yaw) 

CG  -  Center  of  gravity  of  total  weight  of  balloon 


Figure  3. 


Balloon  Geometry  and  Applied  Forces 


equilibrium  angles  in  the  longitudinal  plane  are  not  zero. 
Therefore,  when  solving  the  lateral  degrees  of  freedom,  it 
must  be  assumed  that  the  longitudinal  variables  remain 
constant  and  equal  to  their  equilibrium  values.  In  both  the 
longitudinal  and  lateral  cases,  the  tether  is  simulated  by 
three  rigid  links.  The  number  of  uncoupled  dynamic  equations 
is  four  for  the  longitudinal  response  and  five  for  the  lateral 
response.  j 

B.  GENERAL  STABILITY  THEORY 

The  equilibrium  configuration  of  a  tethered  balloon 
can  be  defined  as  that  position  which  demands  that  the 
summation  of  all  applied  moments  equals  zero.  The  equi¬ 
librium  is  said  to  be  stable  if,  for  any  small  disturbance, 
the  system  ultimately  returns  to  its  equilibrium  conditions. 
Two  types  of  stability  are  of  interest.  In  the  first 
(statically  stable),  a  small  displacement  .f  the  system  wi1! 
create  forces  which  tend  to  return  the  sys  .em  to  its  equi¬ 
librium  position.  The  second  (dynamically  stable)  produces  a 
motion  which  e/entually  restores  equilibrium.  If  the  motion 
is  periodic,  it  is  characterized  by  a  damped  frequency  and  a 
damping  ratio.  Similar  definitions  apply  for  statically 
and  dynamically  ur stable  motions.  A  third  possibility  is  for 
the  system  to  be  neutrally  stable  during  which  the  motion 
neither  diverges  nor  converges. 

Characteristic  equations  were  derived  for  study  of  the 
stability  of  this  system.  The  general  approach  was  as  follows 

(1)  Derive  the  nonlinear  equations  of  motion  in 
three  dimensions  for  each  degree-of-freedom 

(2)  Assume  the  motion  is  near  equilibrium  so  that 
the  aquations  can  be  linearized  and  separated 
into  a  longitudinal  motion  and  a  lateral  motion 

(?)  Laplace  transform  the  linear  equations  from  the 
time  domain  to  the  "S"  domain  assuming  that  the 
initial  conditions  are  zero.  This  establishes  a 
matrix  equation  of  the  following  form: 

l A  j  |  X  (S)(  =  (Of  (1) 

where  {x  (S)}  is  the  eigenvector  and  [A]  is  a  square  matrix 
whose  elements  are  quadratics  in  S  containing  the  physical 
properties  of  the  system. 

(4)  Expand  the  determinant  of  [A)  such  that  the 
characteristic  polynominal  is  obtained. 
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Each  root  of  the  characteristic  equation  represents 
a  term  in  the  general  solution  of  the  form/  A^e  Sit#  where 
is  the  " i"th  root  and  A^  is  an  amplitude,  dependent  on  the 

initial  conditions  of  the  system.  Both  real  and  complex  roots 
may  appear  where  the  complex  roots  occur  in  conjugate  pairs. 

In  general  for  ”n"  degrees  of  freedom,  the  characteristic 
equation  will  yield  "2n"  roots.  Each  pair  of  complex  con¬ 
jugate  roots  represents  one  oscillatory  motion,  while  each  real 
root  represents  one  aperiodic  motion. 


First  consider  an  oscillatory  system.  This  motion  is 
characterized  by  two  roots  of  the  form  S.  =  X  +  i  Y,  where  X 
and  Y  are  real  numbers  and  i  =  J-l.  1  Several  important 
quantities  can  be  found  from  the  root.  The  natural  frequency 

associated  with  this  motion  is  o>n  =  >/x^  +  Y^.  The  damping 

ratio  is  C  =  — ~ —  .  The  damping  frequency  is 

/ - 

w .  =  w  v/  1  ~  C  =  Y.  It  is  also  of  interest  to  know  the  time 
a  n 

to  half  amplitude  for  a  stable  root  or  the  time  to  double 
amplitude  for  an  unstable  root.  This  quantity  can  easily  be 
found  by  considering  one  oscillatory  motion.  The  general 
solution  for  free  vibration  is 

Z  =  Ce  ^  w  n  sin  (ood  t  +  J3)  (2) 

where  0  is  the  phase  angle  dependent  upon  initial  conditions 

C  is  a  constant  dependent  upon  initial  condit  Ions 


The  half  amplitude  time  is 


4.  _  0.693 

21 ' 


(3) 


The  second  possibility  is  an  aperiodic  motion  given 
by  the  expression 

Z  =  CeXt 


where  X  is  the  real  part  of  one  root  and  the  imaginary 
part  (Y)  is  zero 

If  X  is  negative,  Z  approaches  zero  as  time  increases 
indefinitely  and  the  motion  is  said  to  be  overdamped.  Like  the 
oscillatory  motion,  roots  which  give  overdamped  motions  will 
also  occur  in  pairs.  However,  unlike  the  complex  conjugate 
roots  which  result  in  one  oscillatory  motion,  each  real  root 
is  a  distinct  motion.  Therefore,  it  is  possible  for  an  "n" 
degree-of-f reedora  system  to  have  ''2n"  distinct  motions  if  the 
system  is  so  heavily  damped  that  all  the  roots  to  the  charac¬ 
teristic  equation  are  real. 
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There  is  a  third  possible  motion  which  is  a  borderline 

case.  If  two  roots  are  real  and  equal,  the  system  is  said  to 

be  critically  damped.  The  motion  will  be  aperiodic  and  both 
roots  will  give  the  same  motion. 

The  general  solution  to  the  motion  of  the  system  is  a 
linear  combination  of  all  the  motions  defined  by  the  roots  to 
the  characteristic  equation.  Associated  with  each  root  is  a 
mode  shape  which  gives  the  relative  amplitudes  of  each  degree 
of  freedom  when  the  system  is  responding  to  one  particular 
root.  It  is  of  interest  to  establish  these  mode  shapes  so 
that  each  stability  curve  can  be  associated  with  a  definite 

motion  of  the  whole  system.  For  example,  one  mode  shape  may 

show  that  the  pitching  motion  of  the  balloon  is  very  large 
compared  to  the  motion  of  the  tether. 

C.  STABILITY  ANALYSIS 

Derivations  of  the  equations  of  motion  of  the  tethered 
balloon  system  and  development  of  the  characteristic  equations 
for  a  tethered  balloon  system  approximated  with  a  three-link 
tether  are  given  in  Reference  2. 

The  four  linearized  longitudinal  equations  are  Laplace 
transformed,  and  an  eighth  order  characteristic  equation 
generated  which  specifies  stability  characteristics  of  the 
system.  In  like  manner,  the  five  linearized  lateral 
equations  can  be  reduced  to  a  tenth  order  equation  which  gives 
stability  information  in  the  lateral  degrees-of-freedom.  The 
roots  of  these  characteristic  equations  identify  the  natural 
frequencies,  damped  frequencies  and  damping  ratios. 

D.  DYNAMIC  RESPONSE  ANALYSIS 

The  calculations  of  the  balloon  system  response  to 
specific  disturbances  is  the  subject  of  the  dynamic  response 
analysis.  The  most  general  motion  the  system  can  have  is  a 
linear  superposition  of  the  normal  modes. 

Each  aperiodic  or  non-oscillatory  normal  mode  has  one 
arbitrary  constant  {the  initial  value  of  any  one  of  the 
variables)  associated  with  it;  and  each  periodic  or 
oscillatory  normal  mode  has  two  arbitrary  constants  (the 
amplitude  and  phase  angle  of  any  one  of  the  variables) 
associated  with  it.  The  total  number  of  arbitrary  constants 
is  then  eq  nl  to  the  number  of  aperiodic  modes  plus  twice  the 
number  of  periodic  modes;  i.e.  to  the  degree  of  the  character¬ 
istic  equation,  or  the  order  of  the  system.  A  specific 
disturbance  will  excite  the  normal  modes  in  varying  degrees 
and  establish  the  values  of  the  arbitary  constants. 
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The  dynamic  response  of  tethered  balloon  systems  to 
various  wind  disturbances  is  obtained  by  integrating 
numerically  the  longitudinal  and  lateral  equations  of  motion 
to  produce  a  time  history  of  the  dynamics.  The  start 
conditions,  or  equilibrium  conditions  for  the  dynamic  response 
computer  programs  are  obtained  from  the  linearized  stability 
computer  programs  (Reference  2).  This  approach  to  analysis 
has  the  advantage  that  wind  gusts  can  be  produced  and  the 
actual  motion  of  the  system  can  be  observed.  The  major  dis¬ 
advantage  is  that  a  greate”  amount  of  computer  time  is 
required  when  compared  to  evaluation  of  stability  by 
investigating  the  roots  of  the  characteristic  equations. 

The  equations  of  motion  for  the  longitudinal  dynamics 
of  a  tethered  balloon  system,  were  initially  derived  in  two 
forms  (see  Reference  3) : 

(1)  inertia  terms  which  contain  products  of 
angular  velocities  are  neglected, 

(2)  inertia  terms  which  contain  products  of 
angular  velocities  are  included. 

The  concept  of  neglecting  products  of  angular  velocities  is 
associated  with  the  assumption  that  angular  velocities  are 
small;  and  therefore,  products  of  angular  velocities  are 
negligible. 

Numerical  integrations  were  made  with  the  computer  to 
determine  the  effect  of  neglecting  the  inertia  terms  containing 
products  of  angular  velocities.  Although  an  effect  is 
obviously  present,  the  overall  differences  between  the  results 
of  the  two  sets  of  equations  is  small  as  shown  in  Reference  3. 

It  was  decided  that  dynamic  simulation  studies  would 
be  conducted  with  a  model  which  neglects  products  of  angular 
velocities  for  three  reasons.  First,  the  equations  containing 
products  of  angular  velocities  are  shown  to  give  only  slightly 
different  results.  Second,  it  is  desirable  to  keep  the  dynamic 
equations  compatible  to  the  equations  used  in  the  stability 
study  (the  stability  study  used  linearized  equations).  Third, 
it  is  desirable  to  keep  the  longitudinal  equations  compatible 
to  the  lateral  equations  (to  derive  the  lateral  equations  of 
motion  containing  products  of  angular  velocities  in  the 
inertia  terms  would  be  a  very  difficult  task  because  of  the 
number  of  terms  involved. ) 
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SECTION  III 


TEST  BALLOON  DESCRIPTION 


A.  GENERAL 

The  tethered  balloon  system  which  was  flown  to  obtain 
experimental  motion  data  cor. sis ted  of  a  70,000  cubic  foot  aero- 
dynamically  shaped  balloon  produced  by  Lea  Bridge  Industries  of 
Essex,  England,  and  a  Nolaro  tether.  The  Nolaro  tether  was 
0.52  inches  in  diameter  and  weighed  90  pounds  per  1000  feet. 

The  winch  used  is  permanently  installed  at  Fair  Site  in  the 
northwestern  portion  of  the  White  Sands  Missile  Range. 

The  general  arrangement  of  the  70,000  cubic  foot  aero- 
dynamically  shaped  balloon  is  presented  in  Figure  4.  The  major 
elements  of  the  balloon  are  the  envelope,  the  tails  and  the  sus¬ 
pension  system.  The  envelope  is  divided  into  three  separate 
chambers  -  the  qas  chamber,  ballonet,  and  tail  cone.  The  two 
near  horizontal  surfaces  and  the  vertical  rudder  make  up  the 
tails  (Table  I  lists  the  dimensions  for  the  balloon) . 

TABLE  I 

DIMENSIONS  FOR  THE  70,000  FT3  KITE  BALLOON 


PRINCIPAL  PARTICULARS 

Maximum  Capacity  of  Envelopo .  70,000  Ft 

Maximum  Volume  of  Air  Filled  T’cil  (lower  only).  0,000  Ft3 

Maximum  Volume  of  Horizontal  Tails  (each)  .  .  .  3,070  Ft3 

Maximum  Diameter  of  Envelope .  39  Ft 

Length  of  Envelope . 106  Ft  6  In 

Overall  Width  with  Fins  Inflated .  43  Ft 

Overall  Height  from  Bottom  of  Rudder . 48  Ft  6  In 

Height  of  Balloon  from  Rigging  Confluence  Point  66  Ft 

Ballonet  Capacity  .  To  permit  a 

flying  altitude 
of  15,000  Ft 

Bridle  Confluence  Point . 32  Ft  aft  of 


Nose  of  Hull 

45  Ft  below 
Centerline  of 
Hull 
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B.  TEST  BALLOON  GEOMETRY 

The  geometry  of  the  70,000  cubic  foot  balloon 
was  established  with  the  aid  of  data  from  the  balloon 
handbook  (Reference  4)  and  the  data  on  the  BJ  balloon  (Refer¬ 
ence  5) .  The  overall  length  and  diameter  of  the  hull  and 
location  of  the  confluence  point  of  the  suspension  lines  are 
listed  in  Table  I.  The  geometry  of  tails  was  established  on 
the  basis  of  data  in  References  4  and  5. 


The  pertinent  tail  dimensions  for  the  tail  surfaces  are 
the  chord  (A) ,  the  span  (B)  and  the  radius  of  the  tip  (R)  as  de¬ 
fined  on  Figure  4.  From  Reference  5  for  a  BJ  balloon,  the 
proportions  for  the  rudder  (vertical  tail  surface)  relative 
to  the  total  envelope  length  (L)  are 

A  =  0.41  L 
B  =  0.23  L 
R  »  0.135  L 

All  horizontal  fin  surface  dimensions  are  80%  of  rudder  dimen¬ 
sions.  It  is  understood  that  the  tail  surfaces  are  as  used  on 
the  BJ  balloon.  The  test  balloon  has  a  hull  known  as  the  A 
shape  and  this  hull  has  a  fineness  ratio  of  2.7  as  compared  to 
2.5  for  the  BJ  balloon.  Consequently,  for  equal  hull  volume 
the  test  balloon  hull  is  longer  than  the  BJ  balloon  hull.  Tail 
dimensions  for  the  70,000  cubic  foot  barrage  balloon  were  cal¬ 
culated  by  proportioning  the  length  of  an  equivalent  volume  BJ 
balloon.  The  rudder  dimensions  computed  for  the  test  balloon 
are: 

A  =  40  feet 
B  =  22.4  feet 
R  =  13.15  feet 

C.  BALLOON  MASS  CHARACTERISTICS 
1.  General 

The  balloon's  physical  and  apparent  mass  characteristics 
are  presented  in  this  section  and  tabulated  in  Table  II  .  The 
mass  characteristics  include  the  physical  mass,  additional 
masses  in  longitudinal,  lateral  and  vertical  translation,  the 
apparent  masses  in  roll,  pitch  and  yaw;  physical  mass  moments 
of  inertia  and  apparent  mass  moments  of  inertia  in  roll,  pitch, 
and  yaw.  Calculations  include  center  of  gravity,  centers  of 
additional  mass  ana  dynamic  centers  in  pitch  and  yaw. 
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TABLE  II 


BALLOON  MASS  CHARACTERISTICS 


Balloon  Physical  Mass  -  with  payload 

118.10 

slugs 

Center  of  Gravity  -  with  payload  -  Horizontal 

54.44 

feet 

-  Vertical 

6.57 

feet 

Mass  Moment  of  Inertia,  with  payload  -  Roll 

30207 

slug  feet 

-  Pitch 

149022 

slug  feet 

-  Yaw 

134027 

slug  feet 

Product  of  Inertia  -  with  payload 

10391 

slug  feet 

Balloon  Physical  Mass  -  without  payload 

114.04 

slugs 

Center  of  Gravity  -  without  payload  -  Horizontal 

55.24 

feet 

-  Vertical 

5.30 

feet 

Mass  Moment  of  Inertia  -  without  payload  -  Roll 

24913 

slug  feet 

-  Pitch 

141614 

slug  feet 

-  Yaw 

131895 

slug  feet 

Product  of  Inertia  -  without  payload 

7130 

slug  feet 

Additional  Mass  -  Longitudinal 

22.00 

slugs 

-  Lateral 

166.63 

slugs 

-  Vertical 

151.58 

slugs 

Center  of  Additional  Mass,  in  pitch  -  horizontal 

56.18 

feet 

-  vertical 

0.00 

feet 

Center  of  Additional  Mass,  in  yaw  -  horizontal 

58.39 

feet 

-  vertical 

1.66 

feet 

Additional  Mass  Moment  of  Inertia  -  roll 

31422 

slug  feet 

-  pitch 

123811 

slug  feet 

-  yaw 

132268 

slug  feet 

Additional  Product  of  Inertia  -  in  yaw 

-8729 

slug  feet 

Apparent  Mass,  Without  Payload  -  roll 

182.97 

slugs 

-  pitch 

265.62 

slugs 

-  yaw 

280.66 

slugs 

Dynamic  Center,  Without  Payload,  in  pitch  -  horizontal 

55.78 

feet 

-  vertical 

2.28 

feet 

? 

2 

2 

2 


2 

2 

2 

2 


2 

2 

2 

2 


14 


TABLE  IT  (cont.) 
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BALLOON  MASS  CHARACTERISTICS 


Dynamic  Center,  Without  Payload,  in  yaw  -  horizontal  57.11 

-  vertical  3.14 

Apparent  Mass  Moment  of  Inertia,  without  payload 

-  roll  54404 

-  pitch  267256 

-  yaw  264854 

Apparent  Product  of  Inertia,  without  payload  -  in  pitch  7525 

-  in  yaw  -805 

Apparent  Mass,  with  payload  -  pitch  269.69 

-  yaw  284.73 

Dynamic  Center,  with  Payload,  in  Pitch  -  horizontal  55.42 

-  vertical  2.88 

Dynamic  Center,  with  Payload,  in  yaw  -  horizontal  56.75 

-  vertical  3.69 

Apparent  Mass  Moment  of  Inertia,  with  Payload  -  pitch  275832 

-  yaw  267400 

Apparent  Product  of  Inertia,  with  Payload  -  in  pitch  11193 

-  in  yaw  2896 


feet 

feet 

slug  ft2 

slug  ft2 

slug  ft2 

slug  ft2 

slug  ft2 

slugs 

slugs 

feet 

feet 

feet 

feet 

slug  ft2 

slug  ft2 

slug  ft2 
2 

slug  ft 


Notes:  1. 

2. 


Mass  is  in  slugs 

Horizontal  distance  is  feet  aft  of 
theoretical  bow 


3.  Vertical  distance  is  feet  below  hull  £ 

4.  Inertia  is  in  slug  feet2 

5.  Balloon  physical  mass  includes  internal 
gas  and  air 
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Based  on  data  received  from  Air  Force  Cambridge  Research 
Laboratories,  the  estimated  weight  breakdown  for  the  balloon 
system  is  given  in  Table  III.  This  breakdown  is  considered  to 
be  representative  of  the  test  balloon  flown. 

The  payload  includes  batteries,  instrument  package, 
instrument  frame,  cable  cutter  and  load  cell. 

The  density  of  the  air  and  helium  for  the  balloon  mass 
characteristics  analysis  was  based  on  a  temperature  of  50°F. 
Based  on  this  temperature  and  the  reported  cable  tension  (free 
lift)  it  was  estimated  that  the  balloon  was  approximately  88% 
full  of  helium  with  approximately  12%  air  in  the  ballonets. 

2.  Additional  Mass 

The  acceleration  of  the  balloon  in  any  of  the  six  de¬ 
grees  of  freedom  causes  aerodynamic  forces  in  addition  to  the 
velocity  and  attitude  changes.  It  has  been  shown  by  Lamb 
(Reference  6)  that  the  derivatives  of  these  acceleration  forces 
have  the  dimension  of  mass.  For  all  practical  purposes,  this 
additional  mass  term  may  be  added  to  the  actual  mass  of  the 
balloon  for  purposes  of  calculating  forces  and  responses  of  the 
system.  References  7  and  8  were  also  useful  in  the  development 
of  this  section. 

The  additional  mass  and  additional  moments  of  inertia 
for  acceleration  in  a  fluid  have  been  worked  out  theoretically 
for  a  number  of  ellipsoids  of  revolution.  Figures  5  and  6  show 
these  coefficients  plotted  against  fineness  ratio.  The  added 
mass  is  obtained  by  multiplying  these  coefficients  by  the  mass 
of  the  displaced  air.  It  has  been  proposed  and  has  become  the 
custom  to  use  (for  airships  and  balloons)  values  based  on  the 
ellipsoid  having  the  same  volume  and  the  same  length  as  the 
hull.  In  addition,  the  theoretical  longitudinal  coefficient 
of  additional  mass  shown  on  the  curves  are  increased  50%  to 
allow  for  the  boundary  layer  which  is  dragged  along  with  the 
balloon. 

The  fineness  ratio  of  the  equivalent  ellipsoid  to  the 
hull  must  be  computed  based  on  equal  volumes  and  equal  lengths. 

V  =  _IL_  d2  L 
6 

d 


Fineness  ratio 


TABLE  III 


ESTIMATED  BALLOON  WEIGHT  BREAKDOWN 


ITEM 


WEIGHT- POUNDS 


Hull 

Horizontal  Tail 
Vertical  Tail 
Ballonet 
Handling  Lines 
Rip  Panel 

Suspension  Cables  and  Fittings 
Gas  Valves 
Ballonet  Blower 
Pressure  Tubing 
Electric  Cables 
Miscellaneous  Hardware 


Total  Balloon  Weignt 


Payload 


Total  Flight  Weight  (less  tether) 


1543 


1674 


17 


VARIATION  OF  THE  COEFFICIENT  OF  ADDITIONAL  MOMENT  OF 
INERTIA  WITH  FINENESS  RATIO  OF  AN  EQUIVALENT  ELLIPSOID 


Figure  5  .  Variation  of  the  Coefficient  of  Additional  Moment  of 

Inertia  with  Fineness  Ratio  of  an  Equivalent  Ellipsoid 


Figure  6  .  Variation  of  the  Coefficient  of  Additional  Mass  with 
Fineness  Ratio  of  an  Equivalent  Ellipsoid 
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The  flat  tail  surfaces  also  exhibit  an  added  mass  effect 
when  accelerated  normal  to  their  chords.  For  large  aspect 
ratios,  the  added  mass  is  equivalent  to  the  mass  of  air  contained 
in  a  cylinder  whose  diameter  is  equal  to  the  chord  of  the  tail. 
Tests  have  provided  corrections  for  smaller  aspect  ratios  and 
for  taper.  The  tail  added  mass  terms  can  be  obtained  from 
Figure  7  and  added  to  the  hull  added  mass  for  the  total  added 
mass  and  moment  of  inertia  effects,  and  for  determining  the 
dynamic  center. 

The  added  mass  of  the  tail  is  obtained  by  multiplying 
these  coefficients  of  the  mass  of  air  contained  in  a  cylinder 
having  the  length  equal  to  the  span  and  a  diameter  equivalent 
to  the  chord  of  a  rectangular  tail  having  the  same  area. 

3.  Dynamic  Center  Location 

Under  dynamic  loads  the  balloon  will  have  a  mass  and 
moment  of  inertia  which  includes  the  additional  mass  and  addi¬ 
tional  moment  of  inertia  of  the  affected  surrounding  air. 

Because  of  symmetry  in  the  XZ  plane  the  only  product  of 
inertia  considered  is  P  .  Since  the  product  of  inertia  about 
any  axis  of  symmetry  is  zero,  PXY  and  Pyz  are  therefore  zero. 

The  additional  mass  along  each  major  axis  is  added 
directly  to  the  r-ass  of  the  balloon.  The  center  of  the  addi¬ 
tional  mass,  however,  does  not  coincide  with  the  center  of 
gravity.  The  combination  of  these  two  masses  will  determine  the 
dynamic  center  about  which  the  complete  balloon  will  act  in  air. 
Since  the  value  and  location  of  the  additional  mass  is  different 
for  accelerations  along  each  axis,  the  dynamic  center  will  be 
located  at  different  places  for  the  lateral  and  vertical  accel¬ 
erations. 

Figure  8  shows  how  the  dynamic  center  is  located  and  the 
combined  moment  of  inertia  is  calculated. 

D.  AERODYNAMIC  CHARACTERISTICS 

The  static  aerodynamic  characteristics  for  the  test 
balloon  were  developed  from  wind  tunnel  data  obtained  at  the 
University  of  Washington  as  reported  in  References  9,  10  and  11. 
The  aerodynamic  coefficients  as  taken  from  these  references  are 
plotted  in  Figures  9  and  10.  Note  that  the  force  coefficients 
are  based  upon  a  reference  area  of  hull  volume  to  the  two  thirds 
power  which  is  compatible  with  the  equations  of  motion  which 
have  been  developed.  However,  the  moment  coefficients  in  Figures 
9  and  10  are  based  upon  a  reference  length  of  hull  length  and 
are  about  a  reference  center  near  the  bridle  confluence  point  of 
the  bridle  as  noted  in  Figure  10. 
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FIGURE  7.  VARIATION  OF  THE  COEFFICIENTS  OF  ADDITIONAL  MASS  AND  MOMENT 
OF  INERTIA  WITH  ASPECT  RATIO  OF  RECTANGULAR  PLATES 


LOCATION  OF  CA  AND  DC  WILL  BE  DIFFERENT  ABOUT  ALL  THREE  AXES 
M#  -  ADDITIONAL  MASS 
Mg  -  MASS  OF  BALLOON  AND  HELIUM 
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FIGURE  8.  LOCATION  OF  DYNAMIC  CENTER  (DC)  FROM  MASS  CENTER  (MC)  2ND 
CENTER  OF  ADDITIONAL  MASS  (CA) 
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Figure  10.  Lateral  Static  Aerodynamic  Characteristics 


The  static  moment  coefficients  for  the  test  balloon 
have  been  adjusted  to  change  the  reference  length  from  hull 
length  to  and  the  moments  have  been  transferred  from  a 

reference  center  near  the  bridle  confluence  point  to  the  center 
of  volume  of  the  hull.  The  static  aerodynamic  coefficients 
adjusted  for  reference  length  and  reference  center  are  listed 
in  Tables  IV  and  V. 

The  dynamic  aerodynamic  characteristics  of  the  BJ  bal¬ 
loon  were  developed  analytically  as  reported  in  References  1 
and  2.  These  coefficients  were  modified  for  the  test  balloon 
by  accounting  for  the  longer  tail  surface  moement  arms  of 
this  configuration.  The  dynamic  coefficients  are  also  listed 
in  Tables  IV  and  V. 


TABLE  IV 


BARRAGE  BALLOON  LONGITUDINAL  AERODYNAMICS 
(Moments  About  Center  of  Hull  Volume) 

(Reference  Area  V^u^3  "  Reference  Length  V1^3) 


Angle  of  Attack  (a )— 10 

-5 

0 

5 

10 

15 

20 

25 

Lift  Coefficient^. )-.  21 

■L 

-.115 

1 

• 

o 

to 

+.07 

+.16 

+  .255 

+  .345 

+  .47 

Drag  Coefficiently)  .09 

.065 

.05 

.05 

.07 

.11 

.16 

.25 

Pitching  Moment (c  ) — . 280 
in 

-.1397 

-.0352 

+.0595 

+.5194 

-.069 

-.2765 

-.713 

Lift  Force  due  to  Pitch  Velocity  C.  .  <*  2.11/rad 

Pitching  Moment  Due  to  Pitch  Velocity  C  ■ -3. 21/rad 

Drag  Force  due  to  Pitch  Velocity  ■  0.9CTCT  . 


TABLE  V 


BARRAGE  BALLOON  LATERAL  AERODYNAMICS 

(Moments  About  Center  of  Hull  Volume) 

2/3  1/3 

(Reference  Area  "  Reference  Length  V  '  ) 


Static  Aerodynamic 

Coefficients 

Sideslip  Angle  (8) 

-20 

0 

+20 

Side  Force  Coefficient  (C  ) 

+  .450 

o 

• 

o 

-.450 

JT 

Yaw  Moment  Coefficient  (Cn) 

+.0787 

0.0 

-.0787 

Roll  Moment  Coefficient  (C^) 

-.008 

0.0 

+  .008 

Dynamic  Aerodynamic  Coefficients 


Side  Force  Due  to  Yaw  Velocity 

Yaw  Moment  Due  to  Yaw  Velocity 
Roll  Moment  Due  to  Yaw  Velocity 
Side  Force  Due  to  Roll  Velocity 
Yaw  Moment  Due  to  Roll  Velocity 
Roll  Moment  Due  to  Roll  Velocity 


(C  ) 

Yjfc 

*B 

(cn  •  > 
^Bv 

(c..  ) 
,  . 

(cnB’ 

(c"}e' 


+2. 19/Rad 

-3.49/Rad 

-  .442/Rad 
+  .142/Rad 

-  .170/Rad 

-  .327/Rad 


Static  Aerodynamic  Stability  Derivatives 


Side  Force  Due  to  Yaw  Angle 
Side  Force  Due  to  Lateral  Velocity 

Yaw  Moment  Due  to  Yaw  Angle 

Yaw  Moment  Due  to  Lateral  Velocity 

Roll  Moment  Due  to  Yaw  Angle 
Roll  Moment  Due  to  Lateral  Velocity 


(C.  ) 


!4> 


B 


%  1 

(CnvB} 


(C,  ) 


"1 

*B 
<Cju  ) 


B 


+1. 29/Rad 
-1.29/Rad 

+  .226/Rad 

-  .226/Rad 

-  .0229/Rad 


+  .0229/Rad 


SECTION  IV 


FLIGHT  TEST  PROCEDURES 


Flight  tests  were  conducted  to  obtain  experimental  bal¬ 
loon  motion  data  for  comparison  with  analytical  predictions. 

As  stated  previously,  a  70,000  cubic  foot  barrage  balloon  and  a 
0.52  inch  diameter  Nolara  tether  were  flown. 

Instrumentation  consisted  of  the  following  items: 

(a)  Cinetheodolite  coverage  of  balloon  to  obtain 
X,  Y  and  Z  coordinates  of  the  balloon.  Motion 
was  measured  at  the  confluence  point  of 
balloon  suspension  lines. 

(b)  Motion  picture  coverage  of  balloon  motion 
(camera  looked  up  vertically  from  tether 
point  to  give  balloon  yawing  motion) . 

(c)  A  telemetry  package  located  at  the  confluence 
yoint  which  provided 

(1)  Horizontal  relative  wind  speed 

(2)  Vertical  relative  wind  speed 

(3)  X,  Y  and  Z  acceleration 

(4)  Roll  angle 

(5)  Pitch  angle 

(6)  Ambient  pressure 

(7)  Differential  pressure  (helium  compartment) 

(d)  Pilot  balloon  data  which  gives  wind  direction 
and  magnitude. 

Five  flight  test  data  runs  were  made.  The  balloon 
was  displaced  and  then  released  to  obtain  motion  of  the 
tethered  balloon  system.  The  balloon  was  pulled  to  the  side 
to  excite  lateral  motions  and  was  pulled  aft  to  excite 
longitudinal  motion.  A  light  nylon  line  attached  to  the  con¬ 
fluence  point  of  the  balloon  suspension  lines  provided  initial 
displacement  of  the  balloon.  The  line  was  secured  at  the 
ground  a  known  distance  from  the  ground  tether  point.  The 
tether  cable  was  then  payed  out  until  tether  tension  was  re- 
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duced  to  one  half.  A  test  was  initiated  by  cutting  the 
auxiliary  line  at  the  ground. 

The  above  listed  data  items  were  recorded  for  the 
following  test  conditions. 

(a)  Run  #1  -  Lateral  balloon  displacement  of  100 
feet  to  right.  Tether  length  -  1079  feet. 

(b)  Run  #2  -  Longitudinal  balloon  displacement  of  150 
feet  aft.  Tether  length  -  1118  feet. 

(c)  Run  #3  -  Lateral  balloon  fH ^placement  of  150  feet 
to  right.  Tether  length  -  10b 5  feet. 

(d)  Run  #4  -  Lateral  balloon  displc cement  of  150  feet 
to  right.  Tether  length  -  663  feet. 

(e)  Run  #5  -  Longitudinal  balloon  displacement  of  150 
feet  aft.  Tether  length  -  684  feet. 

Table  VI  presents  the  test  log  for  the  flight  test 

program. 
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SECTION  V 


COMPARISON  OF  EXPERIMENTAL  AND  ANALYTICALLY 
PREDICTED  DYNAMIC  BEHAVIOR  OF  TETHERED  BALLOONS 


A.  GENERAL 

The  geometric,  mass  and  aerodynamic  characteristics  of 
the  tethered  balloon  system  which  was  flown  have  been  des¬ 
cribed  earlier  in  this  report.  The  characteristic  equations 
for  longitudinal  and  lateral  motion  were  solved  to  establish 
predicted  natural  frequencies,  damped  frequencies  and  damping 
ratios  for  this  system.  The  dynamic  response  of  the  tethered 
balloon  system  to  disturbances  simulating  test  conditions  was 
obtained  by  numerically  integrating  the  longitudinal  and  the 
lateral  equations  of  Motion. 

Initial  conditions  of  displacement  and  tension  were 
put  into  the  dynamic  simulation  program  as  estimated  from 
flight  test  data  for  a  given  test  case.  During  tests,  an 
auxiliary  load  was  applied  by  the  control  line  at  the  bridle 
confluence  point  to  obtain  initial  balloon  displacement  aft 
or  to  the  right.  For  computations,  the  balloon  was  displaced 
as  in  tests  and  the  reduced  tether  tension  as  a  result  of  the 
auxiliary  load  was  simulated  by  an  artificial  payload  which 
is  removed  at  time  one  second  of  the  calculations. 

Computations  were  made  with  the  stability  and  dynamic 
simulation  computer  programs  to  establish  the  predicted  dynamic 
behavior  of  this  balloon  for  two  test  conditions. 

Longitudinal  Test  Run  No.  5  was  chosen  for  analysis 
because  test  data  indicated  relatively  little  lateral 
motion. 

Test  Run  No.  3  was  chosen  for  the  lateral  motion  test 
since  this  test  had  the  largest  initial  lateral  dis¬ 
placement  and  data  was  available  over  the  greatest 
length  of  time. 

A  comparison  and  correlation  of  analytically  predicted  and  ex¬ 
perimentally  determined  dynamic  response  for  the  tethered 
balloon  system  is  discussed  herein. 

B,  LONGITUDINAL  MOTION 

Calculated  dynamic  behavior  of  the  barrage  balloon  for 
Longitudinal  Test  Run  No.  5  is  plotted  in  Figures  11a,  lib  and 
11c.  An  estimated  wind  speed  of  17.4  feet  per  second  was  ob¬ 
tained  from  the  onboard  cup  anemometer  telemetry  record  (i.e., 
prior  to  release  of  the  balloon) . 
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Figure  11a.  Dyneumic  Response  for  Longitudinal  Test  Run  No. 


or  Longitudinal 


Pilot  balloon  wind  data  was  ratioed  to  this  value  to  obtain  the 
estimated  wind  on  the  tether.  Balloon  pitch  angle  and  fore  and 
aft  motion  as  measured  during  the  tests  are  superimposed  on  these 
plots  for  comparison  purposes.  Correlation  of  pitch  angle  ap¬ 
pears  good.  The  frequency  of  motion  is  comparable  and  damping 
during  the  test  is  somewhat  greater.  Initial  measured  pitching 
motion  may  be  attributable  to  instrumentation  package  motion 
rather  than  the  balloon  itself.  Stability  data  is  listed  in 
Table  VII  for  this  case.  The  first  computed  mode  of  motion  is 
a  damped  mode  which  has  a  damped  frequency  of  .1132  rad/sec. 
Referring  to  the  mode  shapes  of  Table  VII  this  represents  a 
motion  where  the  cable  links  move  as  a  unit  and  balloon  pitch¬ 
ing  motion  is  approximately  180  degrees  out  of  phase.  This 
phase  relationship  is  also  apparent  in  Figure  11a  where  the 
three  links  comprising  the  tether  move  together  and  the  balloon 
pitch  angle  is  180  degrees  out  of  phase.  The  experimental  re¬ 
sults  indicate  a  pitching  frequency  of  0.10  rad/sec  which  is  in 
reasonable  agreement. 

Again  referring  to  Figure  lib,  the  frequency  of  fore 
and  aft  motion  (payload  range)  shows  that  reasonable  agreement 
exists  between  experiment  and  predictions.  Experimental  motion 
is  more  highly  damped.  Although  motion  is  primarily  in  the 
longitudinal  plane,  some  yawing  of  the  balloon  did  occur  during 
the  experiment,  as  shown  in  Figure  12.  This  would  increase 
aerodynamic  drag  and  some  of  the  additional  damping  might  be 
attributed  to  this.  Also,  all  damping  coefficients  used  in  the 
theoretical  predictions  were  calculated  for  zero  angle-of-attack 
and  zero  side  slip  angle.  It  is  apparent  that  the  first  mode  of 
motion  is  the  dominant  one  which  has  been  excited. 

C.  LATERAL  MOTION 

Predicted  dynamic  behavior  of  the  tethered  balloon  sys¬ 
tem  for  Lateral  Test  Run  No.  3  is  plotted  in  Figures  13a,  13b 
and  13c.  Wind  data  used  for  the  computer  simulation  was  based 
on  telemetry  data  from  the  cup  anemometer  before  release  of  the 
balloon.  The  lateral  displacement  (payload  cross  range)  com¬ 
parison  between  prediction  and  experiment  shows  obvious  dis¬ 
crepancies  (Figure  13a) .  Initial  lateral  displacement  is  sim¬ 
ilar  but  as  time  progresses  the  balloon  actually  begins  to  con¬ 
verge  to  a  different  equilibrium  position  and  some  forcing 
function  causes  the  balloon  to  oscillate  with  relatively  high 
amplitude. 

The  stability  analysis  for  Lateral  Test  Run  No.  3  is  listed  in 
Tables  VIII  and  IX.  Referring  to  Table  VIII,  the  first  mode 
of  motion  has  the  following  characteristics: 
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Figure  12.  Measured  Balloon  Yaw  Angle  for  Longitudinal  Test  Run  No 
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.TABLE  IX  (Continued) 


Natural  frequency 


.072  rad /sec 


Damped  frequency 


.069  rad/sec 


Damping  ratio 


.264 


The  modal  analysis {Table  VIII) indicates  that  this  is  a  coupled 
yaw  lateral  displacement  mode  {i.e.,  the  tether  link  angles  a., 
<* 2  an<3  <r  are  approximately  equal) .  Lateral  displacement 
leads  the  balloon  yaw  motion  by  approximately  80  degrees. 

Table  lx  presents  the  longitudinal  stability  characteristics  as 
added  information.  The  experimental  lateral  response  for  Test 
Run  No.  3  is  again  plotted  in  Figure  14  with  additional  data. 
The  frequency  of  this  motion  is  approximately  0.074  rad/sec  as 
noted.  The  ratio  of  successive  amplitudes  of  motion  is 


X1  _  137 
X2  62 


2.21 


The  frequency  of  motion  is  in  relatively  good  agreement  with 
the  small  perturbation  stability  theory  of  .069  rad/sec  damping 
ratio  ( £  )  which  was  predicted  for  the  first  lateral  mode  is 
.264  which  corresponds  to  successive  amplitude  ratio  as 
follows: 


.<5  = 


2  v  £ 

7T~? 


2  x  .264 
v/l  -  (.264)  2 


1.72 


rr~  =  5.59 

*2 

The  analytically  predicted  values  of  amplitude  ratio 
are  substantially  higher  than  observed  in  the  experimental  data 
of  Figure  14.  However,  a  direct  comparison  can  not  be  made 
inasmuch  as  forcing  functions  due  to  wind  gusts  effect  the 
motion.  In  the  stability  theory,  a  free  vibration  condition  is 
assumed. 


A  further  attempt  was  made  to  deduce  the  probable  wind 
fluctuations  with  time  from  the  available  flight  test  data. 

The  relative  wind  which  the  balloon  sees  is  a  vectoral  summa¬ 
tion  of  the  actual  wind  and  the  relative  wind  due  to  balloon 
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Figure  14.  Time  History  of  Lateral  Displacement  for  Test  Run  No 


motion  through  the  air  mass.  Balloon  motion  at  the  suspension 
line  confluence  point  was  observed  by  cinetheodolite.  The 
magnitude  of  the  relative  wind  at  the  confluence  point  was 
measured  by  means  of  a  cup  anemometer  but  direction  of  this 
relative  wi  .id  is  unknown.  Assuming  a  constant  value  of  the 
fore  and  art  wind  relative  to  the  tether  point  and  the  avail¬ 
able  data  an  attempt  was  made  to  calculate  a  possible  lateral 
or  side  gust  time  history  which  might  have  existed  during  the 
test.  A  calculated  balloon  lateral  response  for  this  side 
gust  time  history  did  not  correlate  with  the  experimentally 
determined  displacement. 

In  view  of  the  inability  to  obtain  a  direct  correlation 
between  analytical  predictions  and  experimental  data,  a  more 
fundamental  study  was  undertaken.  This  study  consisted  of  an 
examination  of  the  predicted  response  of  the  tethered  balloon 
system  to  simple  wind  disturbances  and  an  examination  of  the 
effects  of  aerodynamic  damping  on  balloon  motion. 

The  dynamic  response  of  the  balloon  with  a  constant 
wind  gust  of  3  feet  per  second  directed  to  the  right  after  re¬ 
lease  of  the  balloon  is  plotted  in  Figures  15a,  15b  and  15c. 

It  is  apparent  in  Figure  15a  that  the  balloon's  lateral  dis¬ 
placement  is  retarded  and  more  closely  conforms  to  the  actual 
motion  observed  during  the  first  60  seconds  of  the  test.  With 
the  side  gust  maintained  the  balloon  converges  to  a  new  yaw 
angle  and  the  sideslip  angle  approaches  zero  (Figure  15b)  as 
would  be  expected. 

The  nature  of  the  experimental  lateral  displacement 
time  history  for  Lateral  Run  No.  3  changes  at  60  to  90  seconds 
as  observed  in  Figure  14.  It  is  considered  that  a  change  in 
wind  direction  or  magnitude  could  excite  the  lateral  motions 
of  the  tethered  balloon  system.  A  3  fps  side  gust  superimposed 
on  a  steady  state  17  fps  aft  wind  is  approximately  equivalent 
to  a  wind  of  const  ..it  magnitude  but  shifted  in  azimuth  by  10 
degrees.  This  wind  change  could  be  consistent  with  the  light 
and  variable  winds  existing  during  the  time  of  the  test.  Con¬ 
sequently,  it  was  chosen  to  further  investigate  the  dynamic 
behavior  of  the  test  balloon  system. 

The  experimental  yaw  angle  of  the  balloon  during  the 
test  (Figure  16)  can  be  used  to  obtain  initial  conditions  for 
digital  computer  simulation  of  the  mathematical  model.  From 
Figure  16  at  time  63  seconds  yaw  angle  ( )  is  zero  and  the 
yawing  velocity  (  $  )  is  -2.7  degrees  per  second.  (i.e., 
nose  rotating  to  left  when  looking  down  on  the  balloon) .  From 
cinetheodolite  data  the  suspension  line  confluence  point  is 
moving  to  the  right  at  +3.0  fps.  A  dynamic  simulation  of  the 
balloon  motion  was  obtained  with  these  initial  conditions  and 
a  forcing  function  consisting  of  a  side  gust  of  3  fps  to  the 
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Figure  15a.  Dynamic  Response  for  Lateral  Test  Run  No.  3 
Constant  Lateral  Wind  Gust  of  3  fps  to  Right 


FIGURE  16 


SECONDS 


left  with  a  duration  of  40  seconds. 

The  lateral  response  of  the  tethered  balloon  system  to 
this  disturbance  is  duplicated  in  Figures  17a ,  17b  and  17c.  A 
comparison  of  simulated  and  experimental  lateral  balloon  motion 
shows  some  similarities  and  some  differences.  Consider  one  cycle 
of  motion  from  110  seconds  on  after  the  forcing  function  has 
been  removed.  The  simulated  motion  has  a  period  of  98  seconds 
or  a  frequency  of  0.064  rad/sec  and  the  experimental  motion  has 
a  period  of  90  seconds  or  a  frequency  of  0.070  rad/sec.  The.se 
results  are  in  close  agreement.  The  stability  and  modal  analy¬ 
sis  for  this  balloon  system  (Table  IX)  also  identifies  this 
motion  as  the  first  mode  with  a  frequency  of  .069  rad/sec.  The 
mode  is  a  coupled  yaw-lateral  displacement  with  the  lateral  dis¬ 
placement  lagging  the  yawing  motion. 

It  is  also  apparent  in  comparing  experimental  data  and  the  dyn¬ 
amic  simulation  that  the  simulated  motion  at  this  frequency  is 
more  highly  damped.  An  additional  mode  of  motion  is  apparent 
in  the  simulation  after  the  first  mode  has  damped  out.  This 
mode  is  a  slow  convergence  of  balloon  yaw  and  lateral  displace¬ 
ment  to  the  initial  equilibrium  position.  Again,  examination 
of  the  model,  analysis  of  Table  IX  reveals  an  overcritically 
damped  mode  with  balloon  yaw  motion  and  lateral  displacement  in 
phase. 


The  highly  damped  first  mode  indicated  by  the  dynamic 
simulation  is  not  in  accordance  with  the  observed  response. 

This  might  be  attributed  to  the  accuracy  with  which  aerodynamic 
damping  characteristics  can  be  determined  by  calculations  or 
the  inability  to  input  the  proper  forcing  functions.  Consider¬ 
ing  the  former  further  calculations  were  made  with  the  dynamics 
computer  simulation  program  to  determine  the  effects  of  aero¬ 
dynamic  damping  on  system  motion.  The  effect  of  reducing  aero¬ 
dynamic  damping  coefficients  on  lateral  motion  is  shown  in 
Figure  18.  The  lateral  displacement  has  been  obtained  with  a  3 
fps  side  gust  as  before  and  the  six  lateral  aerodynamic  damping 
coefficients  have  been  reduced  to  percentages  of  100%  damping 
presented  in  Section  III.  Although  static  lateral  aerodynamic 
data  were  not  available  for  angles  of  sideslip  greater  than  20 
degrees  a  projection  of  possible  aerodynamic  characteristics 
(Figure  19)  was  made  to  permit  dynamic  simulation  at  larger 
angles.  It  is  apparent  that  lateral  balloon  motion  is  quite 
sensitive  to  aerodynamic  damping. 

In  Figure  18  a  reduction  in  aerodynamic  damping  was 
shown  to  have  a  very  significant  effect  on  the  lateral  dynamics 
of  the  tethered  balloon.  In  fact,  at  a  point  between  81%  and 
85%  of  the  estimated  damping  the  lateral  motion  becomes  unstable. 
Lateral  stability  computer  runs  were  made  with  damping  coeffi¬ 
cients  reduced  by  as  much  as  20%;  and  although  the  roots  of  the 
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ASSUMPTIONS 

1.  LIFT  (CT)  DATA  AVAILABLE  AT  LARGE 
ANGLES  LOF  ATTACK  ARE  INDICATIVE 
OF  NATURE  OF  SIDE  FORCE  (Cy>  . 

2.  YAW  MOMENT  (C  )  AND  ROLL  MOMENT 
(C i)  IS  DIRECTLY  PROPORTIONAL  TO 
SIDE  FORCE  (Cy). 


Figure  19 . 


Estimate  of  Static  Lateral  Aerodynamic  Coeffi¬ 
cients  at  Large  Sideslip  Angles 
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YAW  MOMENT  COEFFICIENT  (C  ) 
ROLL  MOMENT  COEFFICIENT  (C») 


characteristic  equation  become  less  stable,  they  give  no  indica 
tion  of  the  unstable  motion  encountered  in  the  dynamic  simula¬ 
tion.  This  points  out  an  inadequacy  of  the  stability  program. 
Although  a  balloon  can  be  designed  to  be  stable  near  its 
equilibrium  position,  it  may  develop  motions  far  away  from 
equilibrium  which  are  unstable.  It  should  be  pointed  out  that 
the  damping  coefficients  are  considered  to  be  invariant  from 
this  equilibrium  condition.  In  reality  they  would  vary, 
possibly  enough  to  prevent  the  unstable  motion  shown  in  the 
simulation. 


SECTION  VI 


CONCLUSIONS  AND  RECOMMENDATIONS 


Digital  computer  programs  were  developed  previously  to 
describe  tethered  balloon  stability  qualities  and  to  simulate 
dynamic  motions  of  tethered  balloon  systems.  Experimental  test 
conditions  and  test  tethered  balloon  characteristics  were  input 
in  the  computer  programs  and  balloon  dynamic  behavior  was  pre¬ 
dicted.  A  comparison  of  experimental  and  predicted  results 
permits  establishing  the  validity  of  these  mathematical  tools. 
The  following  conclusions  are  drawn  from  this  program. 

1.  Experimental  and  predicted  longitudinal  dynamic  behavior 
are  in  reasonable  agreement.  The  experiment  clearly  excites  the 
first  mode  of  motion  with  a  frequency  as  predicted  by  stability 
theory.  This  mode  is  a  coupled  motion  of  balloon  pitch  and  fore 
and  aft  motion  of  the  tether  as  a  whole  (180  degrees  out  of 
phase) .  Damping  of  the  fore  and  aft  motion  is  greater  in  the 
real  world  and  may  be  attributed  in  part  to  the  fact  that  yaw¬ 
ing  motions  existed  and  would  contribute  to  greater  aerodynamic 
drag  forces. 

2.  The  experiments,  even  though  carefully  controlled  to 
excite  longitudinal  and  lateral  motions  independently,  result 
in  coupled  tethered  balloon  system  motions.  This  may  be  attri¬ 
buted  in  part  to  the  light  and  variable  winds.  It  is  apparent 
that  the  lateral  motion  of  the  test  balloon  in  particular  is 
sensitive  to  wind  changes. 

3.  It  is  more  difficult  to  establish  correlation  of  experi¬ 
mental  and  predicted  lateral  motions  than  longitudinal  motions. 

4.  The  stability  analysis  and  dynamic  simulation  identify 
two  lateral  modes  of  motion  which  are  excited.  The  first 
lateral  mode  is  a  coupled  yaw-lateral  displacement  with  the 
lateral  displacement  lagging  the  yawing  motion.  The  frequency 
of  this  mode  is  in  good  agreement  with  experimental  measure¬ 
ments.  The  second  mode  is  a  highly  damped  mode  with  balloon 
yaw  motion  and  lateral  displacement  in  phase. 

5.  Discrepancies  between  experimental  and  predicted  later¬ 
al  motions  may  be  attributed  to  undefined  wind  forcing  functions 
and  the  accuracy  with  which  aerodynamic  damping  coefficients  can 
be  calculated.  Analysis  indicates  that  lateral  motions  can  be 
relatively  sensitive  to  aerodynamic  damping. 

6.  The  comparison  of  experimental  and  predicted  dynamic 
characteristics  of  the  tethered  balloon  shows  a  reasonable 
match  although  some  discrepancies  do  exist.  In  spite  of  these, 
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stability  and  dynamic  simulation  computer  programs  are  a  power' 
ful  tool  tor  design  and  analysis  of  tethered  balloon  systems. 

As  a  result  of  this  program  the  following  recommenda¬ 
tions  are  made: 

1.  For  future  analysis,  obtain  additional  and  more  accur¬ 
ate  data  for  tethered  balloon  systems  of  interest.  Additional 
aerodynamic  data,  particularly  aerodynamic  damping  data,  may 
be  obtained  with  dynamic  wind  tunnel  models.  Additional 
field  testing  with  accurately  determined  time  histories  of 
wind  magnitude  and  direction  might  also  permit  deduction  of 
aerodynamic  damping  characteristics. 

2.  It  is  recognized  that  some  limitations  exist  in  the 
present  dynamic  simulation  techniques  as  a  result  of  treating 
the  longitudinal  and  lateral  motions  as  uncoupled.  It  is 
recommended  that  the  mathematical  model  be  extended  to  be 
completely  three  dimensional  and  that  provisions  for  a  wind 
vector  changing  in  direction  and  magnitude  as  a  function  of 
time  be  incorporated. 
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